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1.0  INTRODUCTION 


The  work  presented  in  this  report  has  been  directed 
in  two  rather  different  areas,  which  are  discussed  separately 
in  Sections  2.0  and  3.0. 

In  Section  2.0  we  describe  a  new  method  of  computing 
the  dipole  autocorrelation  function  for  molecular  pressure 
broadening.  V7e  believe  that  numerical  implementation  of  this 
approach  should  yield  reliable  estimates  of  far-wing  absorp¬ 
tion  by  H2O  and  other  molecular  species. 

In  Section  3.0  we  document  the  software  which  has 
been  developed  for  updating  the  AFGL  Line  Parameters  Compila¬ 
tion— ^  for  asymmetric  rotor  molecules.  In  Section  3.0  we 
describe  the  function  of  the  various  editing,  sorting,  and 
comparison  subroutines  which  are  used  in  the  updating  proced¬ 
ure.  A  listing  of  the  various  H^O  bands  for  which  the  pro¬ 
cedure  has  already  been  applied  is  provided.  These  updated 
bands  have  all  been  incorporated  into  the  1980  Edition  of  the 
AFGL  Line  Parameters  Compilation. 
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2.0  A  NEW  APPROACH  TO  THE  DIPOLE  AUTOCORRELATION 

FUNCTION  FOR  MOLECULAR  PRESSURE  BROADENING 

The  work  reported  in  this  section  has  benefited 
greatly  through  active  participation  by  S.  A.  Clough  of  Air 
Force  Geophysics  Laboratory,  and  by  R.  H.  Tipping  of  the 
University  of  Nebraska  at  Omaha. 

The  absorption  coefficient  of  a  gas  may  be  computed 
as  the  Fourier  transform  of  the  dipole  autocorrelation 
function—^ 

<j>(t)  =  <y(0)*y(t)>  (1) 

^  iHt  ^  _iHt 

where  y(t)  =  e  h  y  e  is  the  Heisenberg  dipole  moment 
operator,  and  where  the  average  denotes  a  statistical  average 
over  the  canonical  thermal  equilibrium  density  matrix 


p (H)  =  exp(-eH)/Tr{e_BH}.  (2) 

The  absorption  coefficient  can  then  be  expressed  in 
three  equivalent  forms: 

“<“>  •  3SJTnrad.  *"'“>■  (3) 

with 

00 

X"(u)  =  tanh  (^)  /  ||  e"iwt  [<j>Ct>  +  *(-t)],  (4a) 

-00 

.  *  a  -  e-W4")  7  ||  e-lut  *(t),  (4b) 

00 

=  /  ||  e“iwt  [♦(t)  -  ♦c-t)].  (4c) 

*00 


The  equivalence  of  the  above  three  formulas  is  con¬ 
tained  in  the  time-domain  statement  of  the  Fluctuation-Dissi- 

.  3  / 

pation  Theorem-  (FDT), 
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<K-t)  =  <j>(t  +  i  #0, 


(6) 


which  follows  from  cyclic  invariance  in  expression  (1)  for 
the  dipole  correlation  function. 

The  statement  of  the  FDT  in  the  frequency  domain  is 


where 


0°  , 

4>(±w)  =  hr  J  dt  e“lwt  (|>(±t). 


In  near-wing  pressure  broadening  theory,  it  is  not 
terribly  crucial  that  the  FDT  be  satisfied  in  any  approximate 
theory.  This  is  not  the  case  in  considering  far-wing  absorp¬ 
tion,  as  the  following  simple  example  will  illustrate.  In 
particular,  if  one  evaluates  formulas  (4a)  -  (4c)  using  the 
well-known  impact  approximation,—^  one  obtains  the  results 


X"(u>)  =  tanh  (^)  7  p.  |<i  |  |y|  Jf>  |  2 
1  if  x 

.  I  [ - I -  ♦  - I - ], 

*  (ft)-wfi)2+r2  (w+w^)  2+r2 

X n ( w )  =  (1  -  e“Bnw)  y  P.  |<i|  lu I |f> I  2 

if  x 


(u>-wfi)2+r2 


(8a) 


,  (8b) 


X"(w)  =  l  p.  | <i | |v| | f > | 
if  x 


(w-wfi)2+r2  (w+wfi)2+r2 


(8c) 


In  the  extreme  far-wings  (u>  >>  for  all  strong 
lines),  these  three  formulas  predict  very  different  results 
for  the  estimated  absorption.  Thus  we  view  it  as  crucial  that 
any  treatment  of  far-wing  absorption  should  satisfy  the  FDT. 
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In  the  course  of  the  present  contract  we  have  formu¬ 
lated  a  method  of  constructing  the  autocorrelation  function  in 
such  a  way  that  it  satisfies  the  FDT  on  a  microscopic  basis. 
The  approximations  of  the  theory  are: 

a.  The  binary  collision  approximation  (this  is  well-justi¬ 
fied  in  the  far  wings  for  atmospheric  densities). 

b.  The  uncoupled  line  approximation  (the  validity  is  not 
easy  to  assess,  however,  one  might  expect  line-coupling 
effects  to  average  out  in  the  far-wing  region). 

c.  Isotropic  interactions  are  treated  exactly;  anisotropic 
interactions  are  treated  by  second  order  perturbation 
theory.  This  latter  approximation  might  be  improved  on 
by  making  the  approximation  of  ignoring  the  effects  of 
time-ordering  in  treating  various  time-development 
operators. 

A  complete  derivation  of  the  formalism  described 
above  is  presented  in  Appendix  A.  Here  we  briefly  outline  a 
few  of  the  more  relevant  results. 

Within  the  uncoupled  line  and  binary  collision  ap¬ 
proximations,  we  find  that  the  autocorrelation  function  is 
given  by 

4>(t )  =  v  l  p(e.  )  (2  j  .  +  1)  |<j.||y|  |jf>|  2 

hif  1 

•  e  r  (e.  -  e.  )t  C.f(t).  (9) 

*n  it 

This  result  satisfies  the  FDT  provided  C^(-t)  =  (^(t+iMi) . 
For  the  C^(t)  correlation  functions,  we  obtain 

Cif(t)  =  exp  <Np  [q^f(t)  -  1])  (10) 

where  Np  is  the  number  of  perturbers,  and  with 
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=  TTjTTTT  •  (jflm£m| jflj^p 

**>• 

•  Tr  (p(K0)  <3^!  |U(t-iB*r)|jimi> 

•  <jfmf |U(t)+|jfm£>}.  (12) 

In  this  expression,  the  bracketed  quantities  are  Clebsch-Gordon 
coefficients,  and  the  U’s  are  time-development  operators. 

From  Eqs.  (9),  (30),  (11)  one  can  indeed  show  that  the  FDT  is 

ft 

satisfied,  along  with  the  reality  condition  C^(t)  =  C^(-t). 

Going  to  second  order  perturbation  theory  in  the 
anisotropic  interaction,  we  find 

Cqif(t)  -  1]  =  K..(t)  +  Kf(t),  (12) 

with 

Kf(t)  *  -  It  I  »<cj>  o'6!!’ 

•* 


l  l<jfJ^|Vtj*J»^f>|2  f(t). 


(13) 


jiJ’k’ 


where  p(ej),  p(e+)  are  density  matrices  for  the  internal 
states  of  the  perturbers,  and  the  translational  motion  of  the 
relative  coordinate  (respectively),  and  with 


f(t)  *  -  3*-  It  -  1  ~  ) , 

a  6  ap 


(14a) 


Cl-cos(ioagt )  ] 


U) 


ap 


cu 


H  Ccua Bt  •  sinCu>afjt)] 


(14b) 


aB 


with 


_  1 


(t-i  rit  “  l.rl?l  )• 


“aB  =  X  vtjfJE  ”  fcj 


(IS) 
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Similarly 


K  <t)  =  -  I  p(eT)  P <£=►> 

1  *2  jS  J  * 


(16) 


jiJ’k' 


f(t)  =  -i-  { (t  -  iBtt) 
4o6 


a) 

♦  -i-  (e  lwaB 


iw^.e(t-i$}f) 


U) 


-  1)>, 


(17a) 


cxB 


or 


(t)  = 


u 


—  [1  +  B«u>  -  e^aB  cos  w  ftt] 
2  ap  ap 


aB 


+  srr  ^aB1 " eWia,aB  Sin  waBt] 


aB 


waB  E  %  (ejiJlc  “ 


(17b) 

(18) 


For  simplicity,  in  writing  down  the  above  expressions  we  have 
suppressed  Clebsch-Gordon  coefficients  and  m-summations,  which 
are  important  only  for  numerical  calculations. 

For  large  times,  it  may  be  shown  that  the  above 

US/. 

equations  yield  the  impact  approx imat ion-^-t—  (including  both 
width  and  shift),  i.e.  the  argument  of  the  exponential  in  Eq. 
(10)  is  a  linear  function  of  t. 

At  very  small  t,  the  correlation  function  behaves 
essentially  as  a  Gaussian.  This  implies  that  all  time-deriv¬ 
atives  of  C--(t)  exist  at  t  =  0.  This  in  turn  implies  that  all 
0  /  ^ 

moments—  of  the  lineshape  function  (in  the  frequency  domain) 
are  well-defined  and  finite. 

In  order  to  implement  the  above  theory  for  numerical 
computation,  the  following  tasks  need  to  be  performed: 
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a.  The  Clebsch-Gordon  algebra  must  be  done  (this  has  already 
been  carried  out  for  the  case  where  the  anisotropic  in¬ 
teraction  V  is  taken  to  be  a  dipole-dipole  interaction. 

b.  The  relative  coordinate  wave  functions  which 

satisfy 


cTs-  *  =  esi*r*  (19) 

where  VQ  is  the  isotropic  potential,  need  to  be 
calculated. 

c.  Matrix  element  of  the  form 


must  be  evaluated. 


d.  Various  summations  over  quantum  indices  need  to  be  car¬ 
ried  out,  and  a  Fourier  transform  of  Cif  (t)  must  be  per¬ 
formed  to  yield  the  lineshape  function  C^Cw). 

In  previous  reports  we  have  discussed  some  of  the 
problems  associated  with  obtaining  an  accurate  Fourier  trans¬ 
form  of  C^f(t).  In  particular,  one  must  accurately  sum  con¬ 
tributions  of  C^f(t)  from  time  values  which  may  differ  by 
three  orders  of  magnitude  (e.g.  the  difference  between  a  dura¬ 
tion  of  collision  t,,  and  a  time  between  collisions  x  ).  We 

d  c 

have  also  discussed  a  ’’difference  function’*  technique  for  per¬ 
forming  the  Fourier  inversion.  In  this  method,  one  writes 


C(t)  =  CBasis(t) 


+  C 


diff 


(t). 


(20a) 


where 

Cdiff(t)  =  C(t)  "  CBasisCt>*  (20b) 

The  basis  function  is  chosen  to  have  two  basic  properties: 
a.  CBasis^  bas  an  analyt^c  Fourier  transform. 


7 


b*  Cdiff  (t)  vanishes  at  large  t,  eliminating  the  need  to 
integrate  over  large  t  values. 

We  have  previously  reported  application  of  this 
scheme  to  simple  (real)  correlation  functions  such  as  the 
Anderson-Weiss  correlation  function.  Here  we  briefly  present 
some  results  for  a  complex-valued  correlation  function,  with 
a  time-dependent  f(t)  function  given  by  Eq.  (14). 


The  function  we  have  chosen  has  the  form 

f  .t.  S  -  “ClF(t  ) 


C(t )  =  e"a^r;, 
where  a  is  a  constant,  and 


(21) 


F(t )  =  /  /uT  dw,  e~  knT  f(t). 


(22) 


with 


f(t)  = 


l-CQ8(»lk-«lr.ot)t 


sin(uk-w  .)t 


Although  this  model  represents  a  rather  gross  simpli 
fication  of  the  general  theory,  we  note  that  f(t)  has  the  same 
time-dependence  as  in  Eq.  (14),  with  the  identification  uag 


(u»rot  -  w^).  If  we  choose  a  to  have  the  value 
a  =  . r  ■  e^wrot/kBT, 


(24) 


then  T  can  be  interpreted  as  the  impact-approximation  half¬ 
width,  and  the  long-time  behavior  of  C(t)  is  given  approxi¬ 
mately  by 


2 


C(t)  =  e 


■rt 


M.  a) 


rot 


exp{-  ir  t  Ce  ~Jc^T  ( 


V 


1/2 


“rot 

T TT— 


KX1Z. 


-> 


rot 


B  2 

/  dy  ey  ]}.  (25) 

0 


(The  above  expression  ignores  small  terms  in  the  arguments  of 
the  exponentials  which  behave  as  constants  at  large  t.) 


To  apply  the  difference  function  method,  we  have 
modeled  the  basis  function  Cgas^s(t)  such  that  the  real  part 
is  a  convolution  (in  time)  of  a  Gaussian  and  an  exponential, 
and  the  imaginary  part  is  the  derivative  of  s«ch  a  function. 
This  choice  has  the  correct  properties  that  Re  C(t)  is  even, 
while  Im  C(t)  is  odd.  The  adjustable  constants  in  the  con¬ 
volutions  are  chosen  such  that  Cj.£^(t)  ■*  0  for  large  t.  We 
have  also  been  able  to  choose  C^a  .g(t)  such  that  C^iff^*^ 
also  vanishes  at  t  =  0. 


It  turns  out  that  one  can  obtain  an  analytic  approx¬ 
imation  to  the  Fourier  transform  of  Eq.  (21),  which  is  valid 
in  the  far  wing,  by  expanding  to  first  order  in  the  density, 
i.e. 

C(t)  *  1  -  aF(t).  (26) 

j  This  can  be  Fourier-transformed  analytically and  discarding 

spurious  delta  function  contributions  at  the  line  center  w  =  0 
(which  are  improperly  treated  by  the  above  expansion),  one 
obtains  the  result 


n/..\  —  C(ta>) 
R<“>  - 

(0 

=  e-  EJT  ( 


rot 


+w  1/2 


u> 


-) 


rot 


0(wrot+w), 


where  0  is  the  unit  step  function 


(27) 


0(x)  =  1,  x  >  0 
=  0,  x  <  0, 
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and  where 


CLor(w) 


i  r 
*  w2+r2 


is  the  Lorentzian  lineshape  function. 

In  Figures  1  and  2  we  show  plots  of  the  real  and 
imaginary  parts  of  the  difference  function  C^iff^),  as  used 
in  the  numerical  Fourier  transform  of  C(t).  These  functions 
are  plotted  versus  a  length  variable  x  =  2trct,  and  the  main 
point  to  notice  is  that  the  difference  functions  oscillate  to 
zero  at  large  x  (large  t). 

The  results  of  our  numerical  Fourier  transform, 
with  the  choice  of  parameters  =  100  cm”1,  T  =  296°K  are 

shown  in  Figure  3.  To  within  plotting  accuracy,  the  results 
agree  well  with  the  predictions  of  the  analytic  formula  (27). 

The  above  model  is  rather  too  crude  to  assign  much 
physical  significance,  e.g.  in  a  more  realistic  model,  wrot 
would  be  replaced  by  a  whole  distribution  of  frequencies  (both 
positive  and  negative)  determined  by  the  collision  dynamics 
and  multipole  selection  rules.  In  spite  of  these  complica¬ 
tions,  the  results  of  the  model  calculation  provide  some  con¬ 
fidence  that  an  accurate  Fourier  transform  of  realistic  corre¬ 
lation  functions  can  be  obtained  using  the  "difference  func¬ 
tion"  technique. 
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Figure  1.  Real  part  of  the  difference  function 
fHff(t)  plotted  v:s.  the  length  variable  x  = 


m 


I  «CT ( CM ) 


ry  part  of  the  difference 
ion  Cj •££ (t ) • 


Figure  3.  The  ratio  of  the  lineshape  function  C(w)  to  the 
Lorentzian  lineshape  function  for  the  case 
=  100  cm"1,  T  =  296°K. 
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3.0  SOFTWARE  FOR  UPDATING  ASYMMETRIC  ROTORS  ON 


THE  AFGL  LINE  PARAMETERS  ATLAS 


The  procedure  we  have  employed  to  update  HjO  bands 
for  the  AFGL  Line  Parameters  Compilation  is  based  on  the  pro¬ 
cedure  which  we  first  developed  to  treat  the  v2  band  of 
H20(161).  The  development  of  the  method  is  based  on  two  de¬ 
mands;  to  eliminate  entirely  the  making  of  typographical  er¬ 
rors  when  updating  bands  and,  given  an  old  and  a  new  data 
set,  to  construct  from  these  data  the  best  and  most  complete 
data  set  possible  for  addition  to  the  AFGL  line  atlas.  Some 
of  the  problems  that  can  occur  when  dealing  with  very  large 
data  bases  are:  1)  usually  the  new  data  does  not  fully  re¬ 
place  the  old,  thus  the  final  data  set  must  be  a  controlled 
merging  of  the  two  sets;  2)  when  parts  of  the  new  data  set 
are  in  question,  both  new  and  old  data  are  carefully  analyzed 
to  determine  which  is  the  better  and  hence  which  will  be  re¬ 
tained  for  the  next  version  of  the  atlas;  3)  sometimes  mul¬ 
tiple  quantum  assignments  of  the  lines  occur,  when  they  do 
they  must  be  found  and  eliminated.  Once  the  procedure  has 
been  applied  to  a  band,  we  generate  comparisons  of  the  old 
data  to  the  new  data  to  display  the  changes  that  have  occurred. 
We  have  also  found  it  advantageous  to  keep  track  of  the  source 
of  each  parameter  for  both  future  updating  and  for  discussions 
with  the  users  of  the  AFGL  line  atlas. 

Below  we  give  a  general  discussion  of  the  software 
systems  that  make  up  the  procedure,  and  in  Figure  4  we  give  a 
flow  chart  of  the  method.  Although  there  are  some  changes 
that  must  be  made  from  band  to  band  and  from  molecule  to 
molecule,  we  have  attempted  to  keep  the  procedure  as  general 
as  possible  for  asymmetric  rotor  molecules  (i.e.  H20,  0^, 

S02;  etc.).  To  date  the  procedure  has  only  been  applied  to 
vibration-rotation  bands  of  the  water  molecule,  and  in  Table 
1  we  list  the  updated  bands  along  with  pertinent  information 
describing  the  data  sets.  All  the  new  parameters,  i.e.  line 
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Table  1.  Updated  Bands  of  Water 


positions,  intensities,  and  lower  state  energies  of  transi¬ 
tions  for  these  bands  are  taken  from  C.  Camy-Peyret ,  J.  M. 
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Flaud,  et  al,- -  and  we  wish  to  thank  these  authors  for 

making  their  data  available  to  us.  All  of  these  bands  are 
present  on  the  1980  version  of  the  AFGL  Line  Parameters 
Compilation. 

The  start  of  the  procedure  is  dependent  on  the  form 
in  which  the  new  data  is  available  to  us.  If  we  are  not  able 
to  obtain  the  data  directly  from  the  source,  we  construct  hand 
punched  card  files  of  the  published  data.  These  files  must 
be  carefully  checked  for  typographical  errors.  Once  the  new 
file  is  verified  to  match  the  published  results  we  extract 
the  corresponding  band  from  the  AFGL  tape  and  both  files  (new 
and  old)  are  run  through  the  programs  QCHECK  and  IDCHECK 
which  test  quantum  assignments  according  to  selection  rules 
and  band  assignments,  respectively.  From  this  point  on  we 
simply  refer  to  the  data  used  in  updating  the  band  as  the 
"new  data”  and  the  data  from  the  AFGL  tape  as  "old  data." 

Once  we  are  assured  the  data  files  are  correct  we 
frequency-sort  each  file  and  rotational-quantum-number-sort 
(J',  Ka',  Kc',  J",  Ka",  Kc")  each  file.  The  new  quantum- 
number-sorted  file  is  input  to  the  program  MSFGENF,  and  the 
old  quantum-number-sorted  file  is  input  to  the  program  MSFGENA. 
MSFGENF  and  MSFGENA  create  mass  storage  fiZes  of  the  new  and 
old  data  sets  respectively.  These  mass  storage  files  are  con¬ 
structed  by  the  unique  rotational  quantum  number  assignment 
of  each  transition  and  can  reduce  the  effort  of  searching  for 
a  particular  line  by  several  orders  of  magnitude.  The  result¬ 
ing  mass  storage  file  has  as  its  key  the  upper  state  rota¬ 
tional  quantum  numbers,  and  each  record  contains  all  transi¬ 
tions  with  the  same  upper  state  assignments  and  different 
lower  state  assignments.  The  index  key  of  each  record  is 
calculated  by  the  formula  INDX  =  J'(J'+ 1)  +  Ka'  -  Kc’  +  1. 

If  a  file  does  not  contain  any  transitions  belonging  to  a 
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particular  index,  a  flag  is  written  to  that  index  meaning  no 
information  available  for  that  index.  This  check  is  incor¬ 
porated  into  all  programs  which  read  the  mass  storage  files. 

Once  the  mass  storage  files  (MSF's)  are  constructed 
they  are  tested  by  the  program  MULTQ  which  checks  for  multiple 
assignments  which,  due  to  the  structure  of  the  MSF's,  are 
easily  uncovered.  If  any  multiple  assignments  are  found,  the 
extra  lines  are  deleted  from  the  quantum-number- sorted  file, 
the  MSF  is  reconstructed  and  checked  again  by  MULTQ,  and  this 
is  repeated  until  the  final  MSF's  are  known  to  be  correct. 

We  next  work  with  the  new  and  old  frequency  sorted 
files  and  the  new  and  old  mass  storage  files.  As  stated  ear¬ 
lier,  it  is  required  to  have  a  tabulated  listing  of  the  dif¬ 
ferences  between  the  new  and  old  data.  This  is  done  by  the 
program  COMPARE  which  uses  the  frequency-sorted  new  file  and 
the  old  mass  storage  file.  This  yields  a  line-by-line  list¬ 
ing  of  the  difference  in  the  data  sets  along  with  several 
types  of  flags  that  signal  large  differences  in  the  data. 

These  listings  have  proved  very  useful  in  consultations  with 
users  of  the  AFGL  Compilation. 

The  frequency-sorted  new  file  is  also  used  in  the 
program  HWGEN.  The  new  files  for  h^O  usually  contain  the 
transition  frequency,  the  strength  of  the  transition,  the 
lower  state  energy,  and  quantum  identification  of  the  transi¬ 
tion.  Before  the  file  can  be  used  for  the  Atlas  we  need  to 
add  the  air-broadened  half-widths  for  the  transitions.  HWGEN 

adds  to  the  new  file  the  corresponding  half-widths  as  calcu- 

18  / 

lated  by  Davies  and  Oli. —  In  addition,  the  program  also 
adds  ID  labels  to  identify  the  source  of  each  datum  (five  ID 
labels  in  all  -  source  of  the  transition  quantum  numbers,  the 
source  of  the  frequency,  the  strength,  the  half-width,  and 
the  lower  state  energy).  The  ID  becomes  important  in  future 
references  concerning  the  quality  of  the  data  for  a  particular 
line.  The  ID  labeled  file  containing  the  corrected  half¬ 
widths  is  now  referred  to  as  the  frequency-sorted  new  file. 
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The  frequency-sorted  old  file  and  the  new  MSF  are 
used  in  the  program  OLDLINE  which  selects  all  lines  in  the 
old  data  which  are  not  present  in  the  new  data.  This  file, 
which  is  called  the  missing  line  file,  is  then  run  in  the  pro¬ 
gram  TAGMISS.  The  TAGMISS  program  uses  two  MSF’s  containing 
the  rotational  energy  levels  of  the  upper  and  lower  vibra¬ 
tional  states  of  the  band  in  question  to  add  to  the  missing¬ 
line  file  the  line  positions  and  lower  state  energies  derived 
from  the  new  data.  This  is  not  always  possible  since  some  of 
the  energy  levels  for  certain  vibrational  states  are  unknown; 
in  these  cases  the  old  data  is  retained.  The  program  adds  ID 
labels  to  the  lines  according  to  the  types  of  changes  that 
have  occurred.  Finally  all  lines  that  retain  the  old  fre¬ 
quency  are  inspected  and  the  line  position  is  estimated  from 
the  new  data.  This  is  done  to  remove  any  obvious  typograph¬ 
ical  errors  in  the  old  data. 

At  this  stage  of  the  procedure  we  have  an  ID  labeled 
new  file  and  an  ID  labeled  missing  line  file  from  the  old  data. 
These  two  files  are  fed  into  CMPLETE  which,  using  mass  storage 
techniques,  adds  to  the  two  files  whenever  possible  more  accu¬ 
rate  half-widths  computed  by  a  method  similar  to  reference  18 
but  with  a  more  correctly  determined  Anderson  parameter  bQ  for 
the  complete  transition  i  -*•  f.  In  carrying  this  out, 
gram  corrects  the  ID  labels  when  necessary  and  merges  the  two 
files.  The  file  is  frequency  sorted  and  used  in  the  program 
ENDBAND  which  writes  the  final  file  to  output  (replete  with 
ID  labels  for  all  parameters),  and  writes  a  finaj.  file  to 
AFGL  format  for  addition  to  the  next  edition  of  the  AFGL  Line 
Parameter  Atlas.  ENDBAND  also  gives  useful  statistics  on  the 
data  making  up  the  band. 

The  final  file  is  used  with  the  old  MSF  to  generate 
a  series  of  plots  that  visually  display  the  changes  that  re¬ 
sulted  from  the  updating  procedure.  This  presently  includes 
a  plot  of  the  changes  in  the  line  positions  and  several  plots 
of  the  changes  in  the  line  intensities.  Examples  of  the  plots 
are  given  in  Figures  5  to  9. 
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Figure  5.  Comparison  of  1980  v1  file  to  1978  file 

a.  Frequency 

b.  Strong  intensities 
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flFGL  1978 


The  procedure  discussed  above  has  been  successfully 
applied  to  the  14  bands  of  HjO  listed  in  Table  1. 
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The  Dipole  Autocorrelation  Function  for  Molecular 
Pressure  Broadening:  A  Quantum  Theory  Which 
'  Satisfies  the  Fluctuation-Dissipation  Theorem. 
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ABSTRACT 


The  dipole  correlation  function  for  molecular  pressure 
broadening  is  treated  in  a  quantum  theory  which  rigorously 
satisfies  the  Fluctuation  Dissipation  Theorem  on  a  microscopic 
level.  The  two  basic  approximations  in  the  theory  are  the 
uncoupled-line  and  binary  collision  approximations.  In  the 
present  paper,  the  consequences  of  the  formulation  are  anal¬ 
yzed  up  to  second  order  in  the  anisotropic  molecular  interac¬ 
tion.  An  isotropic  potential  is  also  included,  which,  in 
principle,  is  treated  exactly.  At  large  times  it  is  shown 
that  the  theory  reduces  to  the  well-known  impact  approxima¬ 
tion.  At  short  times,  an  autocorrelation  function  of  Gaussian 
form,  with  a  renormalization  of  the  initial  state  occupancy 
is  obtained.  It  is  found  that  the  qualitative  features  dis¬ 
cussed  above  are  unaltered  in  higher  order  perturbation  theory 
The  results  imply  that  all  time-derivatives  of  the  autocorrela 
tion  function  at  t  =  0  exist.  This  further  implies  that  all 
moments  of  the  lineshape  function  in  the  frequency  domain 
exist,  hence  that  the  lineshape  function  must  decay  "exponen¬ 
tially"  sufficiently  far  in  the  wings. 


I.  INTRODUCTION 


In  the  calculation  of  molecular  pressure  broadening,  ex- 

1  3 

tensive  use  has  been  made  of  the  autocorrelation  function 
for  the  dipole  moment  operator  in  the  time  domain.  Using  this 
method,  the  power  spectrum  for  absorbed  or  emitted  radiation 
is  obtained  from  the  Fourier  transform  of  the  autocorrelation 
function.  This  result  follows  from  the  quantum  analog  of  the 
Wiener-Khintchine  theorem1  for  classical  correlation  func- 
tions.  In  spite  of  the  diverse  formulations  which  may 
be  found  in  the  literature,  it  appears  that  certain  aspects  of 
the  theory  remain  incomplete  or  have  not  been  adequately 
discussed. 

It  is  generally  assumed  that  the  far  wings  of  spectral 
transitions  are  strongly  influenced  by  the  small  time  behavior 
of  the  autocorrelation  function.  However,  any  rigorous  treat¬ 
ment  for  small  times  must  certainly  satisfy  the  Fluctuation- 

Dissipation  Theorem  (FDT)  .1’13"’1®  The  FDT  has  been  incorpor- 

17 

ated  into  a  number  of  essentially  phenomenological  theories, 

20  21 

using  the  Egelstaff  *  -Schofield  complex-time  transformation 
These  theories  force  the  FDT  to  be  valid  on  a  line-by-line 
basis,  using  an  autocorrelation  function  involving  a  number  of 
adjustable  parameters,  typically  a  time  between  collisions. 
tc,  and  something  analogous  to  the  duration  of  collision,  x^. 
To  our  knowledge,  no  microscopic  formulation  which  satisfies 
the  FDT  has  been  presented  in  the  literature.  In  this  paper 
we  provide  such  a  formulation.  The  present  theory  will  also 
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provide  an  explicit  microscopic  description  for  the  short-time 
behavior  of  the  autocorrelation  function,  and  for  the  transi¬ 
tion  to  the  long-time  regime  where  the  impact  approximation  is 
generally  considered  to  be  valid.  These  two  aspects  have  re¬ 
mained  unclear  in  most  previous  formulations. 

In  the  present  paper  we  shall  make  essentially  two  approx¬ 
imations;  (a)  the  binary  collision  approximation,  and  (b)  we 

5  fi 

shall  eventually  ignore  line-coupling  effects.  *  For  problems 
involving  atmospheric  densities,  the  binary  collision  approxi¬ 
mation  appears  to  be  justified  from  experimental  observations. 
The  neglect  of  line-coupling  effects  is  more  serious,  however, 
construction  of  a  general  theory  which  includes  such  effects 
rigorously,  and  which  also  satisfies  the  FDT ,  appears  to  be  a 
more  complicated  problem.  Within  the  two  approximations  above, 
we  shall  show,  at  small  times,  that  the  autocorrelation  func¬ 
tion  for  each  line  is  of  Gaussian  form.  This  implies  that  all 
time-derivatives  exist  at  t  =  0,  hence  that  all  moments  of 
the  lineshape  function  must  also  exist.  The  transition  to  the 
long-time  impact  regime  will  also  be  examined. 


i 


34 


II.  GENERAL  THEORY 


form 


We  shall  write  the  absorption  coefficient  (cm”1)  in  the 
17 


4  tt2o)  n 


rad 


a  (to)  =  - -  X"U> 

X” (to)  =  tanh  (-^)  /  e-iu)t  C^(t)  +  4>(-t)] 


(II-l) 

(II-2a) 


=  (1  -  e 


t  dt  -iwt  ,  x 

>  /  7 iff  e  *(t) 


,  (II-2b) 


=  /  ^  e“lwt  C<t>(t)  -  4><-t>]  .  (H-2c) 


In  these  equations,  nrad  is  the  number  density  of  radiating 
molecules,  B  s  (kgT)_l,  and  $(t)  is  the  autocorrelation  func¬ 
tion  given  by 

4>(t)  =  Tr  {p(H)jj(0)»jj(t)}  ,  (II-3a) 

<K-t)  =  Tr  {p(H)y(0)*v(-t)J  ,  (II-3b) 

=  Tr  { p(H)y ( t ) (0 ) }  ,  (II-3c) 


where  p(t)  is  the  Heisenberg  operator 

_  iHt  .  iHt 

w(t)  =  e  M  p(0)  e”  M 


and 


p (H)  =  e"3H/Tr  {e"pn} 

=  .“*H/Z 

is  the  canonical  density  matrix. 


-BH, 


,  (II-4) 

(II-5a) 
,  (II-5b) 
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The  Hamiltonian  H  is  for  a  system  consisting  of  one  radi¬ 
ating  (absorbing)  molecule,  and  Np  perturbers  with  which  it 
may  interact. 

The  equivalence  of  formulas  (II-2a)  -  (II-2c)  is  contained 

17 

in  the  time  domain  statement  of  the  FDT,  i.e. 

<p(-t)  =  4>(t  +  ie«)  .  (II-6) 

This  result  is  readily  proven  from  Eqs.  (II-3a  -  II-3c)  using 
cyclic  invariance  and  the  fact  that  e* iHt AH  connnute  with  p(H); 
Eq.  (II-6)  must  also  be  consistent  with  the  relation 

<J>(-t)  =  <|>(t)*  ,  (II-7 ) 

which  guarantees  that  a (a)  is  real. 

From  Eqs.  (II-3a  -  II-3c)  we  may  write 
iHt  iH  (t-iB#) 

4>(t)  =  Z-1»Tr{y  e  if  *y  e"  if  }  (II-8a) 

♦  C-t)  =  4>(t  +  i&h) 

iH  (t+iWO  .  iHt  . 

=  Z_1«Tr{e~lf  y*e  y)  .  (II-8b) 

In  what  follows  it  will  be  convenient  to  decompose  the 
Hamiltonian  H  as  follows: 

H  =  HQ  +  V  ,  (II-9) 

where  contains  the  unperturbed  energies  of  the  molecules 
and  any  purely  isotropic  interactions,  V^,  which  do  not  in¬ 
volve  the  internal  states,  i.e.  Vq  depends  only  on  center  of 
mass  separations.  The  perturbation  V  is  then  taken  to  contain 
the  anisotropic  interactions  between  the  radiator  and  the  per¬ 
turbing  molecules. 
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Now  in  Eqs.  ( II-8a  -  II-8b)  we  introduce  the  time-develop¬ 


ment  operator  defined  by 

iHt  iH0t 

e“T  =  e”’3?  U(t) 

iHt  ..  lH0t 

e  TR  =  U(t)T  e  *T 

or,  more  generally,  for  complex  t 

iHz  _iHoz 

e  Af  =  e  U(z) 

with  adjoint  relation 


,  (II-10a) 


,  (II-10b) 


,  (II-lla) 


e  Af  -  U(z) 


iH0z 


. (II-llb) 


From  the  above  relations,  one  readily  sees  that  U(t)  is  uni¬ 
tary  for  real  t.  Making  use  of  these  definitions,  we  obtain 


iH0t 


iH0t 


♦  <t)  =  v  Tr{p(Hn)U(t-i&K)y  U(t)  +  *e*3?  y  e~Ti  }  ,(II-12a) 


4>C-t )  =  4>(t  +  i&tt) 


iH0t 


iH0t 


=  v  Tr{p(Hg)  e  -ft  y  e  4?  U(t)»y  U(t-i&h)  }  ,(II-12b) 

where  p(HQ)  =  e-0HO/Tr{e“0HO} ,  is  the  unperturbed  density 
matrix,  and  where  v  is  the  ratio 


_  Z0  _  Tr{e"0HO} 
V  =  '  Trte^} 


.  (11-13) 


Froir.  Eq.  (IT-12b)  it  is  clear,  i-  3  trivial  fashion,  that 
the  FDT  is  still  sati  ;fied  at  this  point.  However,  there  are 
ways  of  rewriting  all  of  these  equations  which  make  identifica- 
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tions  far  less  obvious.  For  example,  if  we  let  t  ■*  -t  in 
Eq.  (IT-12a),  we  obtain 


<j>(-t)  =  v  Tr(p(HQ)  U(-t-iBfi)  y  U(-t)+*e"“E  y  e~X~) .  (11-14) 


iH0t 


iH0t 


The  equivalence  of  Eq.  (II-12b)  and  Eq.  (11-14)  is  now  not 
trivial,  and  it  must  be  obtained  from  various  identities  which 
we  will  presently  derive.  An  even  more  subtle  problem  con¬ 
cerns  the  following  point.  One  does  not  obtain  a  correct  ex¬ 
pression  for  4»(t+i0h)  by  simply  replacing  t  by  (t+i&ft)  in  the 
right-hand  side  of  Eq.  (II-12a).  The  reason  is,  in  the  U(t)^ 
operator  of  Eq.  (II-12a),  the  two  operations  of  complex  time 
translation  t  -*■  t+i&tf  and  of  taking  the  Hermitian  adjoint  do 
not  commute.  We  shall  show  that  the  correct  procedure  is  to 
eliminate  the  adjoint  operation,  and  then  to  perform  the  com¬ 
plex  time  translation. 

The  basic  identity  which  we  shall  need  is  obtained  as 

ff 

follows;  in  Eq.  (II-lla)  we  replace  z  by  -z  and  obtain 

* 


iHz 

e“?T 


iH0z 


=  e  3?  U(-z  ) 

We  now  equate  Lqs.  (II-llb),  (11-15)  and  obtain 

iHoz*  *  iHoz* 

U(z)T  =  e  H  U(-z  )  e"  * 


.  (11-15) 


(11-16) 


Some  useful  special  cases  for  real  t  are 


iH0t 


iH0t 


U(t)T  =  eTT“  U(-t )  e  3T 
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( II-17a ) 


U(-t)T  =  e“*"3t  U(t)  e  -h 


,(II-17b) 


iHOt  x  iHOt 

U(t)  =  eTT”  U(-t)T  e 


. (II-17c) 

The  equivalence  of  Eqs.  (II-12b),  (11-14)  for  can  readily 

be  established  by  starting  from  Eq.  (II-12b),  making  use  of 
Eq.  (11-16)  with  z  =  t  -  i&K,  inserting  Eq.  (II-17c)  for  U(t), 
and  finally  using  cyclic  invariance  of  the  trace  expression. 

Next  we  return  to  the  problem  of  correctly  obtaining 
$(t+iBh)  from  <|>(t).  Starting  from  Eq.  (II-12a),  we  make  use 
of  Eq.  (II-17a)  to  eliminate  the  adjoint  operation.  This  con¬ 
verts  Eq.  (II-12a)  to  read 


iH0t 


iH0t 


<Kt)  =  v  Tr{p(HQ)  U(t-ifJft)  y  e  X  »U(-t)  y  e  M  } 


(11-18) 


The  complex  time  translation  t  t+i&ft  then  give6 


<J>(t+i8h)  =  v  Tr{p (Hq)  U(t)  y 

iHQ(t+i8h)  iHQ(t+ie«) 

•  e - *  U(-t-ifJK)  y  e  *  }  .  (11-19) 

Now  inserting  Eq.  (II-17c)  for  U(t),  Eq.  (11-19)  simplifies  to 

give  our  previous  expression  (11-14)  for  $(t+i0h)  =  <^(-t). 

We  have  discussed  the  above  points  in  some  detail  because 
we  will  need  to  make  use  of  similar  manipulations  tq  verify 
that  the  FDT  is  satisfied  at  later  stages  of  the  formulation. 

For  later  purposes,  it  will  also  prove  convenient  to  in¬ 
troduce  the  integral  equations  for  the  U(t)  operators.  With 
the  boundary  condition  U(0)  =  1,  one  obtains  the  integral 
equations 


, (II-20a) 


U(t)  =1  -  £  /  V(t’)  U(t»)  dt' 


U(t)+  =  1  +  k  /  U(t‘)+  V(t')  dt  * 


,(II-20b) 


with 


iHQt 


iH0t 


V(t )  =  e  'ft  V  e 


More  generally,  for  complex  t,  we  write 


. (II-2Cc) 


U(z)  =  1  -  L  J  V(z'  )  U(z»)  dz1 


,  (II-21a) 


with  V(z’ )  = 


iH0z' 


iH0z' 


, (II-21b) 


where  the  integration  can  be  taken  over  any  path  in  the  com¬ 
plex  z'  plane  where  the  integrand  is  analytic. 

To  conclude  this  section,  it  is  interesting  to  compare 
Eq.  (II-12a)  for  $(t)  with  the  result  one  obtains  in  an  anal¬ 
ogous  theory  which  does  not  satisfy  the  FDT.  In  such  a  theory 
which  ignores  so-called  "back  reaction,”  p(H)  is  approximated 
by  the  unperturbed  density  p(Kg),  and  leads  to  the  result 


iHgt 


iHgt 


♦(t)  *  Tr{p(Hn)  y  U(t)  e“ff  y  e  U(t)} 


.  (11-22) 


Comparison  of  Eqs.  (II-12a),  (11-22)  shows  that  the  more  com¬ 
plete  theory: 


(a)  Contains  the  factor  v 


(b)  Replaces  a  U(t)  factor  by  U(t-i&tf), 


(c)  Contains  a  rearrangement  of  the  various  factors  which' 
cannot  be  eliminated  using  cyclic  invariance  of  the 
trace. 

Items  (a)  and  (b)  above  might  have  been  anticipated,  however, 
item  (c)  is  more  subtle,  and  we  shall  refer  to  it  as  a  "re¬ 
ordering”  effect. 


III.  UNCOUPLED- LINE  AND  BINARY  COLLISION  APPROXIMATIONS 


In  this  section  we  shall  make  the  two  approximations  dis¬ 
cussed  in  the  Introduction,  namely  the  uncoupled-line  and 
binary  collision  approximations.  To  facilitate  these  approxi¬ 
mations,  it  is  convenient  to  separate  out  the  internal  states 
of  the  radiating  molecule  in  the  expression  for  Hq.  Thus  we 
write 


,  (III-l) 


where  H^j  contains  only  the  internal  coordinates  of  the  radi- 
ating  molecule.  Then  Hq  contains  all  the  rest  of  the  unper¬ 
turbed  Hamiltonian,  i.e.  the  internal  coordinates  of  all  per- 
turbers,  the  translational  coordinates  of  all  molecules  (in¬ 
cluding  the  radiator),  and  the  isotropic  interactions  between 
the  radiator  and  perturbers.  The  unperturbed  density  matrix 
then  factors  as  p(HQ)  =  p(Hq).  The  operator  p(Hq)  can 

also  be  further  factored,  however,  it  is  convenient  not  to  use 
this  at  present. 


We  denote  the  eigenstates  of  H^  by  |jm>,  where  n.  is  the 
magnetic  quantum  number,  and  where  j  stands  for  all  other 
quantum  numbers  necessary  to  specify  the  internal  state.  Now 
in  Eq.  (II-12a)  we  take  the  trace  over  p(H^)  and  find 
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(E ♦  r  ~£  *  t  )t 

3f  3i 


*  <3^  |p  I  jim]> 

•  Tr{p(H0Xj  !m!  |U(t-iftK)  |  jin»i>»<jfmf  |U(t)+|  j£m£>MlII-2) 

-Be.,,  r  -Be 

with  p (e .  ,  )  =  3i/J  C2j+1)  e  Bej  .  (III-3) 

3i  3 

In  obtaining  the  above  result  from  Eq.  (II-12a),  we  have  used 


the  fact  that 

***0*  ^  iHRIt  ^  iHRIt 

e  y  e  -K  =  e  5?  y  e~  X.  , 

because  the  remainder  of  the  unperturbed  Hamiltonian  (HQ)  com¬ 
mutes  with  y.  It  is  again  interesting  to  compare  Eq.  (III-2) 
with  the  analogous  result  derived  from  Eq.  (11-22)  which  does 
not  satisfy  the  FDT: 


<J>(t)  = 


l  l 

iijjDfjf  mimimfI! 


»(e .  )  e  £  (e. ,-e. , )t 


•  <jimi|y  |jfmf>*<j^m|.|y|j!m!> 

•  Tr{p(HQ)  <jfmf  |U(t)*|  j£m£>«<j!m!  |U(t)  |  j^iri£>}  .  (III-4) 

Comparison  of  the  above  equations  shows  four  differences: 

(a)  The  factor  of  v  =  Zq/Z  in  Eq.  (III-2), 

(b)  UCt)  U(t-iBH), 

(c)  The  difference  in  two  statistical  factors  p(e.,)  and 

3i 

P  (  E  •  )  j 

Ji 

(d)  The  different  order  of  U  relative  to  11^  in  the  trace 
expressions. 


I 
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Items  (c)  and  (d)  above  arise  from  the  "reordering"  effect 
discussed  previously. 

At  this  point  we  shall  make  the  uncoupled  line  approxima¬ 
tion,  jf  =  jf  and  =  j^.  In  a  separate  publication,  one  of 
the  present  authors  (R.  W.  D.)  will  present  a  somewhat  simpli¬ 
fied  approach  which  is  meant  to  be  applicable  only  in  the  far 
wings.  This  formulation  is  carried  out  in  the  frequency  do¬ 
main,  it  satisfies  the  FDT,  and  it  correctly  includes  all 
possible  line-coupling  effects.  We  note  that  the  uncoupled 
line  approximation  eliminates  part  of  the  "reordering"  effect 

discussed  above,  i.e.  p(e.f)  -*•  p(e.  )  in  this  approximation. 

3i  ]i_ 

If  we  now  apply  the  uncoupled  line  approximation  to  Eq. 

(III-2 ) ,  and  then  make  use  of  the  Wigner-Eckhart  theorem  for 

the  y  matrix  elements,  we  obtain 

<Kt)  =  v  l  p(e  .  )(2j  .+1)  |<5.|  1  p  I  |  jf>|2 
3i3f  1 


e  £  (e.^-e.^t  C.f(t) 


,  ’(III-5 ) 


Jf  Ji 

where  the  reduced  matrix  elements  satisfy  the  symmetry  relation 
(2jf+X)|<jf||M||ji>|2  = 

and  where  C^(t)  is  a  correlation  function  given  by 


nj'Si'v  „  m,l 

rn.mimi.mi 


m 


I  •  111  •  iti  rJll  r-l 

1  1  f  f 

•  ( j  f  lmfm  j  j  f  1  j  )  ( j  f lm£m  |  j  f  1  j  ^mj ) 

•  Tr{p(H0><jim! | U(t-i&h) j j^m^> 

•  <jfmf  |U(t)+|  jfmp) 


.  (III-7) 
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Before  proceeding  to  make  the  binary  collision  approxima¬ 
tion,  it  is  interesting  to  note  the  normalization  properties 
of  Eqs.  (III-5),  (III-7).  First  we  note  that  C^(0)  i  1. 
Rather,  from  Eq.  (III-7)  we  obtain 

Cif(0)  =  (2l T+TT  J,  Tr{p(H0)<jim[|U(-i&K)  |jimi>} 

•  l  (jflmfm| jfljimi)(jflmfm| 5flj£»|)  .  (III-8) 


Performing  the  sum  over  m^,  m  in  Eq.  CIII-8)  simply  produces  a 


factor  of  6  . ,  and  leads  to 


cif(°)  =  CiC0)  =  ^4+1  l  Tr{p(H0)<jimi|U(-iB«)|jimi>}.CIII-9) 


1  m 


Although  this  does  not  equal  unity,  we  can  write 
Cif(t)  =  Ci(0)  cif(t) 


,(111-10) 


C.f(t) 

Cif(t)  =  C±( 0 ) 


,(111-11) 


then  C^f(0)  =  1,  and  we  can  rewrite  Eq.  (III-5)  as 


<J>(t)  =  .  l  P<ej  . )  I  <3  £  I  lul  I  jf>!  2  e  ^  (ej  "ej  . 

D  i  1  f  1  fl 


(HI-12  ) 


where 


p(e .  )  5  v  p(e.  )(2j.  +  l)  C.(0) 
3i  3i  1  1 


. (III-13) 


We  now  show  that  p(eJ  )  is  simply  a  renormalized  initial  state 

3i 

occupancy,  in  particular  that 
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• ( III-14 ) 


To  prove  this,  we  combine  Eqs.  (III-9),  (III-13),  (III-14)  to 
give 

l  p(e .  )  =  v  l  p(e.  )  Tr{p(Hn)<j .m. |U(-i0tf) | j -m.>} 
j.  3i  j.m.  ]i  011  11 

=  v  Tr{p(H0)  U(-iBft) }  ,(111-15) 


I 


where  the  trace  is  now  over  the  complete  unperturbed  density 
matrix.  However,  it  may  readily  be  established  that 

Tr{p(HQ)  U(-iBft)}  =  =  v"1  ,(111-16) 


from  which  the  result  (III-14)  follows  immediately. 

Now  the  binary  collision  approximation  to  Eq.  (III-7)  is 
simply 

NP 

Cif (t)  =  Cqif(t)]  h  ,(111-17) 

where  Np  is  the  number  of  perturbers ,  and  with 

qif(t)  =  rgTTTTJ  I  (jflmfm|  jfljimi)(jflm^n|  jflj^!) 


» 


•  Tr{ps(HQ)<  |Ug(t-i0h)  |  jimi> 


•  <jfmf |Us(t)+| j  > }  . (III-18) 

In  the  expression  (III-18),  the  subscript  S  on  the  p  and  U 

operators  denotes  that  these  operators  now  correspond  to  a 

single  radiator  and  perturber. 

One  next  invokes  the  same  argument  used  by  Baranger^  and 
22  23 

others.  *  Namely,  one  assumes  that  q^*.(t)  has  the  form 


♦ 
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q.f(t)  =  [1  ♦  £  Fif(t>] 


,(III-19a> 


=  Cl  +  Fif(t)] 


,  (III-19b) 


where  fl  is  the  normalization  volume,  and  np  is  the  perturber 
density.  Then  in  the  limit  of  large  Np,  Eqs.  (III-17), 


(III-19 )  yield 


Cif(t)  "  e 


np  Fif(t) 


,  (III-20a) 


C.f(t)  =  e 


Np  Cqif(t)  -  1] 


. (III-?0b) 


In  order  to  justify  the  form  in  Eqs.  (III-19),  we  can 
begin  by  noting  from  analogy  with  Eqs.  (11-20),  (11-21)  that 
the  Ug  operators  in  Eq.  (III-18)  can  be  written  in  the  form 


Us(t)+  =  1  +  ws(t)+ 


U0(t-i&h)  =  1  +  W„(t-ie«) 
o  o 


where 


,  (III-21a) 
,  (III-21b) 


ws(t)  =  -  jj  /  Vg(t')  Ug(t’)  dt» 


,  (III-22a) 


Ws(t)+  =  ~  /  us(t»)+  Vg(t')  dt* 

.  t-i&h 

Wg(t-iWi)  =  -  i  /  Vg(z’)  Ug(z ' )  dz» 


,  (III-22b) 


.  (III-22C) 


Then  in  Eq.  (III-18),  the  term  in  the  product  of  the  U’s  cor¬ 
responding  to  unity  may  readily  be  shown  to  sum  to  unity. 
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i 


which  justifies  the  first  term  in  Eq.  (IJI-19a).  The  volume 
dependence  of  the  second  term  in  Eq.  (III-19a)  will  be  estab¬ 
lished  at  a  later  stage. 

It  is  important,  after  having  made  the  uncoupled  line 
and  binary  collision  approximations,  to  be  able  to  demonstrate 
that  the  FDT  is  still  rigorously  satisfied.  The  proof  that 
this  is,  indeed,  the  case  is  given  in  the  Appendix.  Briefly, 
the  proof  consists  of  first  showing  that  the  FDT  is  satisfied 
provided 

Cf i (-t )  =  Cif(t+i0h)  ,  (III-23) 

(note  the  exchange  of  indices  i,  f).  In  our  binary  collision 
approximation,  (III-20b),  the  condition  (III-23)  is  clearly 
satisfied  provided  we  can  show 

qfi(-t)  =  qif(t+iBJl)  .  (III-24) 

The  details  of  this  proof  are  carried  out  in  the  Appendix.  We 
also  remark,  in  the  uncoupled-line  and  binary  collision  approx¬ 
imations,  that  the  reality  condition,  Eq.  (II-7),  is  satisfied 
provided  C^(t)*  =  C^^(-t),  i.e.  if  q.^(t)*  =  q£f(“t).  This 
relation  can  be  proven  using  manipulations  similar  to  those 
found  in  the  Appendix. 

Returning  to  Eq.  (III-18),  we  can  make  use  of  Eqs.  (HI- 


22)  to  write 


8 


^if^^”^  ~  TTTT+Ty  I  ^  j^lm^rn |  j ^  1  j  .m. )  ( j  - lmpn |  j *1  j  .m! ) 

J  i  numim^m  pn 


•  {Tr{ps(H0)<j.in||Ws(t-ia«)  |j.m.>} 

+  Tr{Pg(H0)<jfmf  |Wg(t)^  |  6m, 


f  »mf 


:m. 

l  l 


+  Tr{pg(H0)<jim!  |Wg(t-i&}\)  I  |Ws<t)+|  jj.mp}}, 

(III-25) 

Once  again,  we  can  compare  this  to  the  result  one  obtains 

from  Eqs.  (11-22),  (III-4)  which  do  not  satisfy  the  FDT.  In 

the  simpler  theory,  Wg(t-i&fi)  is  replaced  by  Wg(t)  in  the 

first  and  third  terms  of  Eq.  (III-25).  In  addition,  in  the 

third  term  of  Eq.  (III-25),  the  order  of  the  initial  and  final 

state  matrix  elements  is  reversed  in  the  simpler  theory.  This 

is  the  only  remaining  "reordering"  effect. 

We  also  remark  that  the  first  two  terms  in  Eq.  (III-25) 

3  24 

roughly  correspond  to  S(b>outer  in  Anderson  theory,  *  while 

the  last  term  corresponds  to  Anderson’s  S(b).  _ .  In  graph- 

inner 

2  5 

ical  perturbation  theory  the  two  types  of  terms  correspond 
to  self-energy  and  vertex  corrections,  respectively. 


49 


IV.  INVESTIGATION  OF  TIME-DEPENDENCE  IN  LOW-ORDER 
PERTURBATION  THEORY 

It  is  clear,  since  the  Ug  operators  in  Eqs.  (III-22) 
satisfy  integral  equations,  that  it  will  be  extremely  diffi¬ 
cult  to  rigorously  treat  the  anisotropic  interaction  V  to  an 
arbitrary  order  in  perturbation  theory.  In  this  section  we 
will  therefore  confine  our  discussion  of  the  theory  to  second 
order.  However,  as  we  shall  point  out,  only  the  quantitative 
details  are  altered  if  the  theory  is  carried  to  third  or 
I  higher  order.  Moreover,  we  also  point  out  that  the  present 

theory  satisfies  the  FDT  order-by-order  in  perturbation  theory 
In  this  section,  in  order  to  provide  a  more  simplified 
|  account,  we  shall  treat  the  internal  states  of  both  the  radi¬ 

ator  and  perturber  as  non-degenerate,  i.e.  we  shall  ignore  the 
(2j+l)-fold  degeneracy  in  the  m  indices.  This  also  means  that 
{  we  shall  suppress  the  Clebsch-Gordan  coefficients  and  the  m 

summations  in  Eqs.  (III-18),  (III-25),  and  shall  set  factors 
of  (2j+l)  equal  to  unity.  These  details  can  easily  be  re- 

<  tained,  but  are  important  only  for  numerical  calculations, 
which  will  not  be  attempted  in  this  paper.  As  a  final  simpli¬ 
fication  in  notation,  we  shall  also  drop  the  subscript  S  on 

<  the  Pg,  Ug,  Wg,  Vg  operators  which  appear  in  Eqs.  CIII-18)- 
<111-25 ),  understanding  that  these  operators  now  refer  to  a 
single  radiator  and  perturber. 

*  We  shall  begin  by  considering  the  second  term  in  Eq.  (IIT 

25),  which  we  now  simplify  to  read 
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(IV-1) 


Kf(t)  =  Tr  (p(H0)  <jf |W(t)+[jf>} 

This  term  enters  in  an  identical  fashion  in  both  the  theory 
which  satisfies  the  FDT ,  and  the  one  which  does  not.  However, 
before  considering  this  term  in  detail,  we  briefly  remark  on 
the  other  two  terms  of  Eq.  (III-25),  which  we  also  simplify 
as 

KA(t)  =  Tr  { p(HQ )  <ji|W(t-iBW)|ji>}  ,  (IV-2) 

K.f(t)  =  Tr  {p(H0)<ji|W(t-ie«) | ji>*<jf |W(t)+| jf>}  .  <IV>3 ) 

Concerning  Eq.  (IV-2),  at  large  times  (|t|  >>  Bfi)  this  term 
can  be  approximated  by 

K.(t)  =  Tr { p (Hq )  <ji|W(t)|ji>)  ,  (IV-4) 

and  the  discussion  in  this  regime  is  essentially  identical  to 
that  which  we  present  for  Kf(t)  at  large  t.  At  very  small 
times,  K^(t)  approaches  the  finite  constant 

K±C0)  =  Tr(p(H0)  <ji|W(-iWi)|fi>}  ,  (IV-5) 

and,  as  discussed  in  Section  III,  leads  to  a  renormalization 
of  the  initial  state  occupancy.  The  discussion  of  Eq.  ( IV— 3 ) 
is  somewhat  more  complicated ,  and  this  will  be  given 
separately. 

In  order  to  perform  the  trace  in  Eq.  (IV-1),  we  need  +-he 
resolution  of  the  identity  operator  in  terms  of  the  product 
states  of  Hg.  This  can  be  written 


I 


IJ’k' 


=  1 


(IV-6) 
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* 


(IV-7) 


With  =  1  I J f>  I  V’j**  > 

where  |jjr>  denotes  an  internal  state  of  the  radiator,  |J’> 
denotes  an  internal  state  of  the  perturber,  and  where 
is  an  eigenstate  of  the  Hamiltonian 


H 


ISO 


*  Vr> 


(IV-8 ) 


In  Eq.  (IV-8),  V2  =  V^,  with  r  =  r^-i^  t^ie  relative  coordin- 

r 

ate,  m  =  n^n^/ (m^+n^ )  the  reduced  mass,  and  VQ(r)  is  the  iso¬ 
tropic  pair  potential.  Since  y,  Vq,  V  do  not  depend  on  the 
center  of  mass  coordinate,  5,  of  the  pair,  it  can  be  seen 
that  the  translational  motion  of  the  center  of  mass  plays  no 
role  in  the  subsequent  theory.  We  thus  have 

•  (IV-9) 

For  unbound  states  1 i|j£>  can  be  taken  as  any  suitably  normal¬ 
ized  set  of  continuum  eigenstates,  and  e£  =  -ft2k2/2m.  However, 
if  VQ(r)  leads  to  bound  states,  sums  over  |i|»£>  implicitly  con¬ 
tain  a  sum  over  the  bound  states.  The  unperturbed  energy 
associated  with  a  state  can  then  be  written 


e 


j 


.  (IV-10) 


If  we  now  evaluate  Eq.  (IV-1)  in  first  order  perturbation 
theory  using  Eq.  (III-22b),  we  obtain 

41)(t>  =  jp  I  P(ej)-<jfJ^|V|jfJ^>  ,  (IV-11) 

icJ 


with 
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.  (IV-12) 


p(c£>  =  e“Be£/J  e“BG£" 

icM 

The  expression  (IV-11)  is  seen  to  be  a  linear  function  of  t 
without  any  approximations.  We  note,  however,  that  K^^(t)  is 
pure  imaginary,  and  therefore  it  contributes  a  phase  shift, 
rather  than  damping,  to  the  autocorrelation  function  C££.(t). 

It  can  also  be  shown,  for  multipole  interactions,  that  K^^(t) 
vanishes  when  a  sum  over  magnetic  quantum  numbers  is  per¬ 
formed2^  (which  is  suppressed  here).  In  what  follows,  we 
shall  therefore  neglect  these  first  order  phase  shifts. 

We  next  evaluate  Kf(t)  in  second  order  perturbation 
theory,  and  find 

K<2)(t)  =  -  ij-  I  pCej)  p(e£)-  I  |  <a  |  V  |  0>  |  2  f(t)  ,  (IV-13) 

where,  for  simplicity,  |a>  =  I  P>  -  an<* 

where 

i(ea-Eg)t 

«*>  ■  *  -  *  Ce'^e  •  aV-1U> 

a  8  a  p 

An  alternative  and  useful  form  for  the  f(t)  function  is 


II 

— —  [1  -  cos(w_at)] 

2  OP 

uo6 

- 

• 

1  !  Kb*  ‘  sln^waB*^ 

,  (IV-15 ) 

“OB 

where , 

a:  e  =  (e,  •Efi)/fi.  In  this  seco.id  fc-rr 
ap  a  P 

it  is  immediately 

clear  that  the  real  part  of  f(t)  is  an  even  function  of  t, 
while  the  imaginary  part  is  odd.  Also  from  Eqs.  (IV-13), 
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(IV-15),  we  see  that  the  real  part  of  '(t)  is  always  nega¬ 
tive  (gives  damping). 

The  time  dependence  of  (t)  for  long  and  short  times 
can  now  be  deduced  from  a  study  of  the  f(t)  function.  First 
we  note  that  f(t)  is  a  perfectly  well-behaved  function  when 
(ea-Eg)  *►  0.  By  expanding  the  exponential  in  Eq.  (IV-14), 
we  find 


f(t)  =  t2  as  -*■  0 


.  (IV-16 ) 


Furthermore  it  is  clear  that  the  above  result  is  pre¬ 
cisely  what  we  get  in  the  small  time  limit  t  ■*  0.  Hence,  for 
short  times, 

K<2><t>  *  -  I  oUj)  p(cE).  I 

^  Jlc  -a  *  .T  *  * 

d  3fJ  K  (IV-17 ) 


Thus  as  t  ♦  0,  Kv  ;(t)  is  a  real,  quadratic  function  of  time, 
and  from  Eq.  (III-20b),  it  leads  to  a  correlation  function  of 
Gaussian  form. 

In  order  to  examine  the  long  time  behavior,  we  make  use 
of  the  identities2*5 


tim  [ 1-cos (w  «t)] 

|t|~  - -  •  »  «<«a-)|t| 

waB 


,  (IV-18) 


|t|>»  -i—  [«  et-sin(w  pt)3  *  t  &- 
woB  aB 


.  (IV-19) 


Hence  the  long  time  limit  of  '(t)  is  given  by 


i 
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K^25(t)  =  -  l  l  p(ej)p(e£) 

Jic 

•  l<jfJ*£|V|j£J*^,>|2  *  £,) 

~  IT  ^  E  p(ej}  p(ej<) 

Jk  j^J'k' 


Pr 


.  <IV-2b) 


•jfJk  "j^J’k' 

The  above  result  is  identical  to  the  impact  approxima- 
tion  *  in  second  order  perturbation  theory.  The  appearance 
of  J 1 1  in  the  damping  term  is  also  familiar  in  the  impact 
theory.  If  one  assumes  Eq.  (IV-20)  to  be  valid  for  all  times, 
including  the  neighborhood  of  t  =  0,  the  factor  |t|  leads  to 
singularities  in  the  derivatives  of  C^f(t)  at  t  *  0,  This 
immediately  implies  that  the  higher  order  moments  of  the 
lineshape  function  (the  Fourier  transform  of  C^Ct))  do  not 
exist  in  the  impact  theory.  As  we  have  seen  from  Eq.  (IV-17), 
the  correct  damping  of  C^f(t)  at  small  t  is  Gaussian,  which 
implies  that  all  moments  of  the  lineshape  function  are  well- 
defined  and  finite.  In  the  frequency  domain,  these  results 
imply  that  the  extreme  far  wings  must  decay  in  some  "exponen¬ 
tial"  fashion,  rather  than  the  simple  (w-ok^)”2  decay  predicted 
by  the  impact  theory. 

For  intermediate  times,  it  is  clear  that  the  time  depend- 
ence  of  K^.  (t)  is  complicated,  and  it  will  probably  need  to 

be  extracted  numerically.  A  rough  criterion  for  the  inter- 


mediate  time  region  is  <u»ag>  t  »  1,  where  is  some  aver- 


dB 


age  energy  difference.  For  a  situation  in  which  *<wag>sJcBT* 
we  estimate  the  transition  region  to  occur  at  times  of  order 
t  *  -tf/kgT  =  Bft  =  2.6  x  lo”1*  sec,  for  T  =  296K.  It  is  then 
clear  that  this  is  precisely  the  range  where  the  difference 
between  W(t-iBfi)  and  W(t)  begins  to  be  important  in  the  K^(t) 
term. 


Before  going  on,  we  remark  that  the  long  time  behavior 
O') 

of  the  term  (t)  (with  W(t-iBft)  a  W(t))  can  be  obtained 
from  Eq.  (IV-20)  with  the  substitutions  -*•  j^,  j^,  and 

with  the  imaginary  term  in  Eq.  (IV-20)  changing  sign.  This 
implies  that  the  real  parts  of  Kp^(t)  and  K^2*(t)  add,  while 
the  imaginary  parts  subtract.  This  is  also  familiar  from  the 
impact  theory. 

A  second  item  that  can  be  disposed  of  at  this  point  is 
justification  of  the  normalization  volume  dependence  in  the 

/o\ 

second  term  of  Eq.  (III-19a),  i.e.  we  now  show  that  (t), 

as  given  by  Eqs.  (IV-17),  (IV-20),  is  correctly  proportional 
to  1/ft.  To  see  this  we  rewrite  Eq.  (IV-12)  as 


p(e£>  =  £  e'Be£/(£  l  e‘6e£") 


.  (IV-21) 


In  the  limit  of  large  ft,  i  7  — - —  /  d*ic",  and  the  remain- 

0  p,  (2ir )  * 

ing  factor  of  ^  in  Eq.  (IV-21)  turns  out  to  give  the  expected 
1  (2 ) 

jy  volume  dependence  of  (t).  To  complete  the  proof,  we 
note  that  if  |i|»£>  is  a  continuum  wavef unction,  it  will  contain 
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a  normalization  factor  of  ,  while  bound  state  |^1>  are  nor- 

/JT  5 

malized  independent  of  fl.  Then  in  Eqs.  (IV-17),  (IV-20),  if 
and  |'l’j£»>  are  both  continuum  states,  the  square  of  the 
matrix  element  yields  a  factor  1/fi*.  This  is  just  what  is 
needed  to  turn  the  Jc,  1c *  sums  into  integrals.  Likewise  one 
sees  that  the  volume  dependence  is  correct  for  bound  ■+•  free, 
and  bound  bound  matrix  elements. 

We  next  consider  the  more  complicated  quantity  lC^(t)  as 
given  by  Eq.  (IV-3).  It  is  also  interesting  to  compare  this 
with  a  quantity 

Kif(t)  =  Tr{p(H0Xjf  |WCt)+|  jf>*<ji|W(t)|ji>}  ,  (IV-22) 

which  obtains  in  the  analogous  theory  where  the  FDT  is  not 
satisfied.  In  comparing  Eqs.  C/-3),  (IV-22)  we  note  that  a 
"reordering"  effect  remains.  However,  a  careful  examination 
of  these  equations  (including  magnetic  quantum  numbers  sup¬ 
pressed  here)  shows  that,  to  second  order,  the  "reordering" 
effect  plays  no  role,  and  only  the  difference  W(t)  W(t-iBfi) 
is  important.  We  also  remark  that  if  one  confines  the  analy¬ 
sis  to  second  order  perturbation  theory,  there  are  a  number  of 
interesting  cases  where  and  vanish.  This  happens  for 
linear  and  asymmetric  rotor  (radiating)  molecules  for  the 
case  of  dipole-dipole  or  dipole-quadrupole  interactions,  be¬ 
cause  one  obtains  diagonal  reduced  matrix  elements  of  the 
radiator’s  dipole  moment  operator.  Inese  results  are  also 
familiar  from  Anderson  theory,  *  i.e.  s(b)inner  makes  no 
contribution  in  such  cases. 


1 
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We  note  from  Eqs.  (IV-3),  (IV-22)  that  Kif(t)  and  Kif(t) 
are  already  at  least  of  second  order  in  V.  Then  to  second 
order,  evaluation  of  the  formulas  for  large  t  gives 

K^Ct)  *  K^Ct)  =  l  l  p(£j)  p(e£) 

Jic  J'£' 


.  (ZV-23) 


This  is,  again,  the  impact  result,  and  (t)  subtracts  from 

(2)  (2) 

the  damping  given  by  the  K)  (t)  and  (t)  terms. 

Next  for  short  times,  Eq.  (IV-22)  evaluates  to  give 

K$f)(t)  =  1  l  P(et)  p(e£> 


Jk  J*k' 


.  (IV-24) 


This  term  is  quadratic  in  time  and  subtracts  from  the  damping 
given  by  (t).  However,  at  very  short  times  we  should 
consider,  instead  of  K^(t),  the  function  K^(t)  which  satis¬ 
fies  the  FDT.  Its  small  time  limit  is 

K^)(t)  =  l  l  p(e  )  p(£-0  .  [e&KUJi^"t0J,J,)  ♦  13 
2*a  jk  j’p 

-ztK  K  p(eJ}  9U& 

41  Jk  J'k* 

.  -  -JL—  - 1] 

•jr-vs* 

•  <3iJ^|V|jiJ,^f><jfJ,^<  |V|jfJ^>  .  CIV-25) 
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I 


1 

! 

1 


i 


i 

s 


I 


I 


(2  ) 

We  note,  although  K>f'(t)  contains  a  term  linear  in  t,  that 

this  term  corresponds  to  a  phase  shift,  rather  than  damping. 

Hence  we  again  find  that  the  damping  occurs  quadrat ically  at 

small  t.  Finally,  for  completeness,  we  quote  the  general  re- 
(2 )  (25 

suits  for  K^£.  (t)  and  K)  (t)  in  second  order  perturbation 
theory; 

Kif}(t)  =  -  l  l  p ( e j ) p ( e £ ) •  < j ± ( V | j £ J ' , > < j f J * i|> £ ,  1 V | j f J 

Jk  J'k’ 

,  Ce  J  k  Jk  -  1] 


[eA<eje~Cj'g')(t'i6K).i] 


.  (IV-26) 


The  long  time  (|t|  >>  WO  limit  (Eq.  (IV-23)  is  obtained  from 
this  formula  by  again  using  the  identity  (IV-18).  The  short 
time  limit  leads  directly  to  Eq.  (IV-25).  For  K^^(t),  we 
find 

K<2)(t)  =  -  ^  l  p(ej)p(e£)  l  l^i^IVlj f(t) 

^  .Tic  -i  r  .1 f  ic  * 

Jk  3iJ  k  ,  (IV-27 ) 

with 


i  (e  -e„) (t-i&H) 

?(t)  ‘  {<t■i6,5,  + 1,5 — tv*? - — 


,  (IV-28 ) 


and  where  eQ  =  j£,  Eg  =  ej*j'jc?*  It  is  not  difficult  to 

explicitly  show  that  our  second-order  result  for  Cq^Ct)  -  1], 
as  obtained  from  Eqs.  (IV-13),  (IV-27),  (IV-26),  satisfies 
the  FDT,  and  it  appears  that  the  FDT  is  satisfied  order-by¬ 
order  in  perturbation  theory. 
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V.  DISCUSSION 


Within  the  binary  collision  and  uncoup led- line  approxima¬ 
tions,  we  have  derived  an  expression  for  the  dipole  auto¬ 
correlation  function  which  rigorously  satisfies  the  Fluctua¬ 
tion  Dissipation  Theorem.  We  have  shown  how  these  results  go 
into  the  impact  theory  at  large  times.  For  short  times  we 
have  found  that  C^(t)  is  damped  in  a  Gaussian  fashion.  Al¬ 
though  our  analysis  in  Section  IV  was  confined  to  second  order 
perturbation  theory  in  the  anisotropic  interaction,  by  going 
to  third  order  we  can  easily  convince  ourselves  that  the  above 
behavior  is  completely  general,  and  only  the  details  are  al¬ 
tered  in  higher  order  perturbation  theory.  These  results  im¬ 
ply  that  the  lineshape  function  in  the  frequency  domain  must 
decay  in  some  "exponential"  fashion  sufficiently  far  in  the 
wings. 

Although  we  have  refrained  from  labeling  the  present 
formulation  as  a  "unified"  theory,  we  believe  that  it  quali¬ 
fies  as  such.  One  reason  that  we  have  avoided  the  above  label 
is  that  most  authors  of  "unified"  treatments  are  content  to 
show  that  their  formulation  goes  into  the  impact  limit  at 
large  times,  and  for  small  times  takes  the  form  of  the  quasi- 

oo.in 

static  or  statistical  theory. 

We  have  shown  that  the  present  theory,  indeed,  does  re¬ 
duce  to  the  impact  approximation  at  large  times.  However,  it 
is  not  clear  that  the  small  tine  limit  of  the  present  formula- 
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tion  has  very  much  in  common  with  the  standard  statistical 
theory. 

One  problem  in  establishing  such  a  connection  is  that 

most  formulations  of  the  statistical  theory  are  what  Smith, 

31 

et  al  have  termed  "scalar"  theories;  i.e.  at  some  conven¬ 
ient  stage,  the  dependence  of  the  potential  on  the  internal 
states  is  ignored,  except  possibly  for  a  constant  which  may 
depend  on  the  various  vibrational  or  electronic  bands  of  in¬ 
terest.  In  our  present  formulation,  if  the  anisotropic  inter¬ 
action  vanishes,  C.^(t)  -  1,  and  there  is  neither  broadening 
nor  shift  of  the  spectral  transition.  This  result  has  also 
been  proven  by  Baranger5  in  a  slightly  different  context, 
within  the  impact  approximation. 

Another  aspect  of  the  simple  statistical  theory  also 
deserves  comment.  In  a  very  terse,  but  highly  illuminating 
paper,  Yakimets3^  has  applied  a  "scalar"  statistical  theory 
to  far  wing  pressure  broadening.  By  invoking  the  high  energy 
approximation,  and  using  the  method  of  stationary  phase, 
Yakimets  derives  a  result  for  the  far  wing  lineshape  function 
for  potentials  of  the  form  (o/rn).  Similar  results  have  been 
obtained  by  Holstein3**  and  others.  The  result  for  the  (pseudo) 
dipole-dipole  case  (n  «  3)  is  that  the  far  wing  decays  as 
(«-«££)“*,  precisely  the  same  dependence  as  the  impact  approx¬ 
imation  would  give.  Although  the  above  result  may  have  valid¬ 
ity  somewhere  in  the  wings,  our  present  formulation  shows  that 
it  cannot  possibly  be  correct  in  the  extreme  far  wings,  which 


must  be  "exponential.”  It  is  clear  that  theorists  such  as 
Yakimets  and  Holstein  are  well  aware  of  the  above  limitations 
however,  there  appears  to  be  a  misconception  among  some 
workers  that  the  statistical  theory  is  applicable  in  the  ex¬ 
treme  wing  region. 

Finally,  although  our  formulation  has  been  completely 
quantum  mechanical,  this  does  not  appear  to  preclude  the  use 
of  semiclassical  methods.  In  particular,  the  eigenstates 
|i|>£>  of  the  isotropic  Hamiltonian  (IV-8)  might  be  chosen  to 
be  time-independent  WKB  wavefunctions .  We  hope  to  explore 
this  and  other  computational  possibilities  in  a  future 
publication. 


APPENDIX 


PROOF  THAT  THE  FDT  IS  SATISFIED  IN  THE  BINARY  COLLISION 
AND  UNCOUPLED- LINE  APPROXIMATIONS 


The  first  step  in  the  proof  is  to  show,  for  an  uncoupled- 
line  approximation  of  the  form  Eq.  (III-5),  that  a  sufficient 
condition  for  the  FDT  to  be  satisfied  is 

Cf.(-t)  =  Cif(t+i0K)  .  (A-l) 

The  above  relation  applies  to  any  uncoupled-line  approximation, 
and  is  not  limited  to  the  binary  collision  case. 

Now,  starting  from  Eq.  (III-5),  we  have 

<K-t)  =  v  l  p(c.  )(2j.+l)|<j.|  |p|  I  jf>|2 

-i  A  J-}  1  1  1 

3i3f  1 

i  (£j  ~£j  )  t 

e"JT  3f  3i  •  (A- 2 ) 

We  next  change  names  of  dummy  variables  i  ♦  f,  and  make  use  of 
the  symmetry  relation  (III-6)  for  the  reduced  matrix  elements. 
This  gives, 


♦  C-t)  =  v  l  pCe.  )(23.+l)|<ji| Jul |jf>|* 
ji3f  f 


i  (e.  -e.  )t 
•  e  3f  3i  Cf.(-t) 


(A-3) 


Retaining  this  result,  we  next  obtain,  from  Eq.  (III-5), 
the  result, 
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♦  (t+iwo  =  v  l  p(e .  )(2j.+l)|<j.|!y|Jj  > 

kif  1 

Cif(t+i&«) 


i  <e.  -e .  )(t+i&ft) 
e#  3f  3i 


(A-4) 


From  Eq.  (III-3),  we  note  that 


p(e.  )  e 
3i 

This  gives 
♦(t+iBh)  =  v 


-8(e.  -e.  ) 


(A- 5) 


I.  p(e^  )(2ji+l)|<ji| lv| |jf>|2 
'i3f  f 


i  (e.  -e.  )t 

e  IT  3f  3i  Cif(t+iBfi) 


(A- 6) 


Comparing  Eqs.  (A-3),  (A-6 )  we  see  that  a  sufficient  condition 
that  the  FDT  theorem  (Eq.  (II-6))  be  satisfied  is  just  Eq. 
(A-l) .  It  is  then  clear,  in  our  binary  collision  approxima¬ 
tion  (III-20b),  that  (A-l)  is  satisfied  provided  we  can  show 
qfi (-t )  =  qif(t+iBh)  .  (A-7) 

To  prove  the  above  result,  we  start  from  Eq.  (III-18) 
and  construct  qf^(-t).  To  carry  this  out,  we  let  t  ♦  -t, 

♦  j^,  and  it  is  also  convenient  to  make  the  following 
changes  in  the  dummy  m-summation  indices: 


m. 


ra 


mj 


m. 


f  -i 
m!  m., 


m 


i  ^  mf 
m  ♦  -m 


(A-8) 
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This  leads  to  the  result 
1 


W-0  *  T77T  (jil»!-m|jiljfBp<j1H.i-»|ji15f»f) 


m^mun^mpn 


•  Tr{ps(H0)<jfmf |Us(-t-ieft)|jfm^> 

•  <jimilUs(-t)+| .  (A-9) 

24  35 

We  next  make  use  of  the  symmetry  properties  *  of  Clebsch- 
Gordan  coefficients  to  obtain  the  identity 
<  j  ilm ! -m  |  j  il  j  fmp  ( j  ilmi-m  1  j  ) 

(2jf+l) 

*  'C73  . V IT  •  (A-10) 

Upon  inserting  this  result  into  Eq.  (A-9),  we  obtain 

i  mumim-mim 

•  Tr{ps(HQ)<jfmf  |Ug(-t-ig«)  Ijfinp 

•  <j£m||Ug(-t)+| j^m^>)  .  (A-ll) 

Retaining  this  result,  we  now  want  to  compute  q^Ct+iWl), 
starting  from  Eq.  (III-18)  for  q^Ct).  Recalling  our  discus¬ 
sion  in  Section  II,  we  have  to  be  careful  in  performing  this 
exercise.  In  particular,  we  have  to  eliminate  the  adjoint 
operation  in  Eq.  (III-18),  before  making  the  complex  time 
translation  t  t  +  i&h.  Thus  we  make  use  of  Eq.  (II-17a) 
to  write 
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-cjfiiif  |Us(t)+|jfm»> 


<jfjnf  I 


iH0t 


Us(-t)  e' 


iH0t 


iH0* 


iHQt 


=  e  -tt  <jfmf |Ug(-t) | jf®f>  e“  H  .  CA-12) 

Now  inserting  (A-12)  into  Eq.  (III-18),  and  then  letting 
t  -*■  t+i&ft,  we  find 

•  Tr{ps(H0Xjim!|Us(t)|jimi> 

•* 

iH0 

•  e  W  (t+i(Jfi)  <3f®f |Ug(~t-i&H) | jf“f> 


iHo 

•  e  3T”  (t+iMO) 

Application  of  cyclic  invariance  gives 


(A-13) 


Tr{  }  =  Tr{e‘ 


iH0t 


<jim||us(t)| ji,ni>  e 


iH0t 


•  Pg(H0)<jfmf  |UgC-t-ieh)  | 


(A-14) 


We  then  note  that 


iH0t  i»0t 

e  H  <ji»i|Us(t) | jimi>  e  M 

iH0t 

s  <j.m!|e  E  Uc(t)  e  "Ti  |j.m.> 

1  X  O  1  X 

s  <3^1  |Us<-t)+|jimi>  ,  (A-15) 

where  we  have  made  use  of  Eq.  (II-17b)  in  the  last  step. 
Inserting  (A-15)  into  (A-14)  yields 
Tr  {  }  = 

Tr  {<jim!^|Us(-t)+| jimi>  pg(H0)<jfmf |Ug(-t-iBft) J jjmp) 

=  Tr  {ps<H0)<jfmf.|Ug(-t-ieft)  j  jfm^> 

•  <jimi|Us(-t)+|jimi>}  .  (A-16) 

Upon  inserting  (A-16)  into  (A-13),  we  find 

qif(t+iBtt)  =  I  (jfrlmfinl  jfljimi)(jflm»m|  jflj.m[) 

Ji  m . ro Im^m^m 

•  Tr{pg(H(J)<jfmf  |Ug(-t-iP«)  |  jjin^ 

•  <j  ^m||Ug(-t)+|  j  £»»£> )  •  (A-17) 

The  above  result  is  seen  to  be  identical  to  Eq.  (A-ll)  for 
qfi(-t).  This  completes  the  proof. 
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